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Abstract: We present a precise theoretical prediction for the signal-background interfer¬ 
ence process of gg{-^ h*) —)• ZZ, which is useful to constrain the Higgs boson decay width 
and to measure Higgs couplings to the SM particles. The approximate NNLO A-factor 
is in the range of 2.05 — 2.45 (1.85 — 2.25), depending on Mzz, at the 8 (13) TeV LHC. 
And the soft gluon resummation can increase the approximate NNLO result by about 10% 
at both the 8 TeV and 13 TeV LHC. The theoretical uncertainties including the scale, 
uncalculated multi-loop amplitudes of the background and PDF-l-as are roughly 0(10%) 
at NNLL'. We also confirm that the approximate iL-factors in the interference and the 
pure signal processes are the same. 
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1 Introduction 


A scalar particle of a mass about 125 GeV compatible with the standard model (SM) Higgs 
boson has been discovered recently at the LHC [1, 2]. It is necessary to pin down its various 
quantum numbers and couplings in order to determine its identity. The total decay width 
of the Higgs boson is an important variable that would appear in all the global fitting 
procedures. If it is measured to be the value predicted by the SM, then the confidence to 
consider this particle as the SM Higgs boson is increased. On the other hand, if it is larger 
than the value predicted by the SM, there would be new decay channels for this Higgs 
boson, e.g., invisible decay channels, or the couplings between the Higgs boson with SM 
particles should be modified. A precise measurement of the total width may open another 
window on new physics. 

However, given that the width of the Higgs boson (~ 4 MeV) is much smaller than the 
energy resolution of the detector (~ 1 GeV) , it is impossible to precisely measure the line 
shape and thus the total width of the Higgs boson at a hadron collider. And one can not 
obtain the total decay width from global fitting of various on-shell production and decay 
channels [3]. Taking the golden channel gg ^ h ^ ZZ as an example, the cross section 
can be expressed as 


9ggh9hZZ 

ruhTh 


( 1 . 1 ) 


where gggh and ghzz denote the couplings between the Higgs boson and other SM particles, 
and Tfi is the total decay width of the Higgs boson. Here we have used the narrow width 
approximation for the on-shell Higgs boson production and decay. It is obvions that a 
simultaneous rescaling of couplings and width would result in the same cross section, which 
means that there is no way to get independent information on the couplings or the width 
from only these kinds of measnrements. 
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Figure 1. The leading order Feynman diagrams for the signal and background. Diagrams (a) and 
(6) denote the amplitude of signal and background, respectively. Their interference contributes to 
the signal. Diagrams (c) represents the dominant SM background. 


Other proposals have been presented in the literatures so far to bound the total width, 
which depend on additional model or mass resolution assumptions [3-8]. For example, 
it is assumed that the Higgs coupling to a IH or Z boson pair is not much larger than 
in the SM [9], or the Higgs couples only to the SM particles [10]. A method to bound 
the Higgs boson width through apparent mass shift due to interference between the Higgs 
resonance in gluon fusion and the continuum background amplitude for gg' —)• /i —)• 77 
was investigated [4, 5], but the experimental mass resolution is modeled by a Gaussian 
distribution for simplicity [ 6 ]. A recent proposal makes use of the information of the cross 
section at the non-resonance region where the final states have an invariant mass larger than 
125 GeV [3, 7, 8 ]. Actually, looking into the non-resonance region, where s ~ 

{Mzz is the invariant mass of the Z boson pair), brings two changes to the formula in 
eq.(l.l). First, the narrow width approximation is not applicable any more and thus 
in the propagator can be neglected. And because of the Z boson pair and top quark pair 
threshold effects, the cross section away from the Higgs threshold is not negligible small. 
Second, the interference process between the signal (figure 1(a)) and background (figure 
1 (b)) becomes important. And the cross section at the non-resonance region is given by 

2 2 

dcT ^ 9ggh9hZZ _ 9ggh9hZZ 

d^zz i^zz ~ i^zz ~ ''^D^zz 

where the first and second term arise from the pure signal and signal-background inter¬ 
ference processes, respectively. From eq.(1.2), we see that the cross section in the non¬ 
resonance region is only sensitive to the Higgs boson couplings. Gombining the information 
from both the on-shell and off-shell regions provides a way to measure or bound the Higgs 
boson total width. 

The above statement is actually also based on assumptions. For example, the couplings 
gggh and ghzz do not vary when changing from on-shell to off-shell regions, or have a 
known dependence on Mzz- And the discussion of the limitations of the off-shell coupling 
measurements has been available recently [11-13]. Nevertheless, with this method, the 
ATLAS and CMS collaborations have obtained the upper limit of F/i < 4.8 ~ T.TF^’^ and 
F/i < 5.4r™, respectively, at a 95% confidence level [14, 15]. 

Moreover, the off-shell Higgs production and decay has significant impact on search 
for new physics in addition to the interpretation of the Higgs total width [16-18]. And 
the on-shell Higgs production can not distinguish the contributions from the htt and hgg 


- 2 - 














(induced by new colored particle loop) couplings since they would give rise to the same 
effective operator for a single Higgs boson on-shell production production and decay. The 
off-shell Higgs production breaks this degeneracy because the hgg coupling is sensitive to 
the off-shellness of the Higgs boson Therefore, from the off-shell Higgs production, one 
can obtain a constraint on the htt coupling with the similar accuracy to that from the 
pp —>• tth production [19]. 

Though the theoretical prediction for the pure signal process is known up to next- 
to-leading order (NLO) and next-to-next-to-leading logarithmic (NNLL) level with the 
finite top quark mass [20-22], the interference process is predicted only at leading order 
(LO)[7, 8, 23]. Some recent papers provide the QCD radiative corrections to the production 
of 55 —)• ZZ via a top quark loop in the heavy top quark limit [24] and via a massless quark 
loop [25, 26]. A high-order soft-collinear approximation has been obtained for the signal- 
background interference of a heavy Higgs boson (600 GeV) production and decay into WW 
in ref. [27], which shows that the A'-factors of the approximate NLO and next-to-next-to- 
leading order (NNLO) corrections are very sizable, and the theoretical uncertainties are still 
large, about 9% at approximate NNLO. This result implies that the soft gluon effects are 
important, which should be resummed up to all order, and provide more reliable theoretical 
predictions. 

There are several differences between a heavy Higgs boson (600 GeV) production with 
decay into WW and the SM Higgs boson (125 GeV) off-shell production with decay into 
ZZ. First, the SM Higgs boson is relative light, and the contribution from the signal- 
background interference is different, compared to the heavy Higgs boson. More explicitly, 
the interference cross section for a heavy Higgs boson is positive while the one for the SM 
Higgs boson is negative. Second, the interference for a heavy Higgs boson is most significant 
around its mass threshold, while the interference for the SM Higgs boson dominates around 
Mzz = 200 GeV. Third, there are additional contributions at high orders in the SM Higgs 
boson (125 GeV) off-shell production and decay into ZZ because of the contributions from 
the ggZ triangle loop diagrams; see the diagram 2(c). Last, searching for WW final states 
requires an additional jet veto applied to suppress the large tt background. This would 
induce another kind of large logarithms that need also to be resummed to all order [28]. 

Note that in the large invariant mass region, the Z bosons are significantly boosted. 
The dominant contribution comes from the longitude component of the Z boson, which is 
similar to a Higgs boson. Therefore, it is expected the impact from soft gluon resummation 
in this process is similar to that in the Higgs boson pair production, which we have studied 
earlier [29]. 

This paper is organized as follows. In section 2, we describe the resummation formalism 
in this process briefly. We then investigate the NLO and NNLO expansions from the 
resummation formalism in section 3. In section 4, we discuss the numerical results including 
the invariant mass distributions and the theoretical uncertainties. The conclusions are given 
in section 5. 

^Another method to break the degeneracy is boosted Higgs production. 
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Figure 2. The next-to-leading order Feynman diagrams for the interference process. Diagrams 
(a) and (6) denote the virtual correction while diagrams (e) and (/) denote the real correction. 
Diagram (c) represents a new contribution, compared with gg —>■ WW, to the background that 
could interference with the signal. 


2 Resummation formalism 


In this paper, we concentrate on the interference process, i.e. the interference between 
diagrams 1(a) and 1(6), which contributes to the signal. The amplitude squared of diagram 
1(a) is taken to be the pure signal, which has been calculated very precisely. The amplitude 
squared of the diagram 1(6) is considered as background. However, this kind of background 
is much less than the one coming from diagram 1(c). The amplitudes of the diagrams in 
figure 1 have been computed at LO, and can be found in gg2VV [7, 23] and MCFM [8]. 

More precise predictions require the calculation of the diagrams in figure 2. This has 
not been achieved so far due to the complex two-loop diagrams with massive particles in 
both loops and external states. Given that the invariant mass of the final state is large, 
the soft gluon contribution is expected to dominate the higher order corrections. And this 
special contribution can be resummed to all orders in Og- In this section, we describe the 
neccessary resummation formalism briefly. A more detailed discussion of the factorization 
and resummation formalism can be found in the single Higgs [30] or the double Higgs [29] 
productions at the LHC. 

In the off-shell Higgs production and decay to an on-shell Z boson pair, the invariant 
mass of the Z boson pair Mzz is so large that additional real emissions are strongly con¬ 
strained. The imbalance between the virtual and real corrections induces large logarithms 
of the form ln”'(l — z), where the partonic threshold variable z is defined as 


z = 



( 2 . 1 ) 


with the variable y/§ being the partonic center-of-mass energy, 
region is defined as 


Mly 

T = 1 


The hadronic threshold 
( 2 . 2 ) 
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with y/s being the collider energy. In the threshold region, the differential cross section 
can be written as 

where ffgg is the luminosity of the initial-state gluons in the protons, defined as 

dx 

ffggiy^fj')= —fg/NAx,^^)fg/N2iy/x^^^)■ (2-4) 

Jy X 

And the partonic differential cross section is given by 

daB{Mzz) = ^'\^d^2. (2.5) 

Here, |A4|^ is the LO color and spin sum (averaged) amplitude squared, and d ^2 denotes 
the two-body phase space of the Z boson pair. In addition, C{z, Mzz, y-f) represents 
the hard scattering kernel. After applying the similar derivation procedure for Drell-Yan 
process [31] in the framework of soft-collinear effective theory (SCET) [32-36], it could be 
factorized to two parts, i.e., the hard function and the soft function, 

C{z,Mzz,yf) = 'HiMzz,yf)S{'/I{l - z),yf). (2.6) 

The hard functions for the signal and interference processes are given respectively by 

TdsigiMzz, y) = \Ch{iMzz,y)\^, 

TdintiMzz, y) = Tie. [Ch{i Mzz, y)Cvv{i Mzz,y)], (2.7) 


where the hard Wilson coefficient Ca{iMzzi y) is obtained by matching the gluon operator 
from QCD to SCET [30]. Here the index a = h and VV accounts for the process of 
gg ^ h* ^ ZZ and gg —)■ ZZ, respectively. The renormalization group (RG) equation for 
the hard Wilson coefficient is 


dCa{iMzz, y) 
din /i 


• cusp 


(a.) In 

y 


+ 7^®(as) 


CaiiMzz, y), 


( 2 . 8 ) 


and the corresponding solution (evolved to the factorization scale) is 


Ca{iMzz,yf) = exp 


—A/p 

‘^S{yh,yf) - ar(//h,Ii/)ln-^ - psa{yh,y}) 

yf 


Ca{iMzZ, yh), 
(2.9) 


where the intrinsic hard scale yq is chosen as y^ = s-iid the vr^-enhanced terms 

are resummed to all order in a^ by RG evolution from = —M^z to positive values of 
y'j [30]. S{u,y) and are functions respectively of the anomalous dimensions T^gp 

and 7 , as defined in [36], and we do not write them explicitly here. At the fixed orders, 
the hard Wilson coefficients have perturbative expressions in series of the strong coupling, 

CX> 

Cai/Mzz,//) = 1 + Y.cirHL) (^)", ( 2 . 10 ) 

n=l 
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where L = In (—In the matching procedure for gg ^ h* ^ ZZ, the Wilson 
coefficient at NLO can be expressed as 


= - CaL'^ + 


TT 


( 2 . 11 ) 


Here we have written the scale dependent part explicitly. ci is scale independent and 
its precise value requires the calculation of the two-loop diagrams, such as the diagram 
2(a), and can be extracted from refs. [37-40]. In the large top quark mass limit, ci = 11. 
In our numerical discussion, we have included the exact finite quark mass effect in the 
signal amplitude. Up to NNLO, i.e. three-loop level, we can only make use of the result 
obtained in the large top quark mass limit at the moment, and cj^ (L) can be approximately 
expressed as 


cf (L) =C\Ha{L) + CATpnfHfiL) + 


TT 


+ Cl 


Hi{L)+C2, 


where 


C2 = 


/7451 2177r2 vr^ 


■ 


54 


499 . 

12 ”8 


and Ha, Hf and Hi have the form 


11 rs .2 f 67 7r2\ /386 llvr^ 

Ha{L)- — + -L +L - 2C^ 

Hj{L) = --L^ + -L^ + L + -vr^ 

9 9 V 279 

Hi{L) = —CaL^ + L ^ 


( 2 . 12 ) 


(2.13) 


(2.14) 


respectively. At the moment, there are no complete results for the two-loop or three-loop 
virtual corrections to the background gg ^ ZZ amplitude. However we solve the RG 
group for the hard Wilson coefficient, and get the exact scale dependent terms as 

41(1) =-CaL^ + 5i, 

c%{L) = CIHa{L) + CATpnfHfiL) + 5iHi{L) + 82 , (2.15) 


where di and 82 represent the unknown scale independent terms at the two-loop and three- 
loop levels respectively. 

The soft function describes soft interactions between all external colored particles. It 
has the same form as in the single Higgs production [30]. Up to 0{as), it is given by 


S (^'\/I(l — z),g^ = <5(1 


^) + 




TT 



5(1 — z) + 6 D{z) 


(2.16) 


with 


D{z) = 


1 


1 - z 




(2.17) 
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This kind of plus distribution comes from the integration with the pole subtracted. The 
soft function obeys an integro-differential evolution equation, which has been shown in 
ref. [29]. Using the Laplace transformation [30], we can obtain the corresponding solution 
as 

S — z), fifj = Us{f^si IJ^f) s{dri, Us) ^ 

where the auxiliary parameter rj is defined as rj = 2ar(Ais,/i/), and 

= exp [-AS{^is,^if) + 2a^w{fis,IJ‘f)] (2.19) 

with ^s being the intrinsic soft scale and usually set numerically. We choose it according to 
the method in [29]. The function s{dri,^s) is the Laplace transformed soft function, which 
is defined as 

/•OO 

s(L,/i,)= / ( 2 . 20 ) 

Jo 

Using the property of RG invariance of the total cross section, we could get the anomalous 
dimension as 

= ^ + j99 + 2j^. ( 2 . 21 ) 

as 

Combining the above components together, we obtain the resummed hard scattering 
coefficient for the interference process 


C{z, Mzz,IJ^f) ='R-e [Chii Mzz, fJ-f)CYv{'i Mzz,fJ-f)] U{Mzz, /^/) 

'(S \ ^-2'rr] 


- M|,(l-.)^ „ ^e 


(1 - z)l- 2 ^ 




+ drji fJis 


r( 2 r?)’ 


( 2 . 22 ) 


where 


aiif^f) 


_\ ~2ar(/J'/i,/J's) 

Z Z 


X |exp [45(/i/j,//5) - 2a^99(/i,j,/is) + 4a.^B(/is,/i/)] I . (2.23) 


This is the main formula in our calculation. In practice, we would use the three-loop 
cusp anomalous dimension and two-loop normal anomalous dimension, and thus denote 
the precision of the resummed result as NNLL', where the prime means that the results 
are not exact NNLL due to the existence of unknown scale independent terms in the hard 
function. Notice that the vr^-enhanced terms have been resummed to NNLL order. Since 
the fixed-order result is only exactly known to LO, we do not match the resummed result 
with fixed-order ones. 
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3 NLO and NNLO expansions 


In the resummed cross section, there are three scales, fif, fj-s and If we set them equal 
to each other, then we obtain the threshold singular contributions, which should appear in 
the hxed-order calculations. Up to NNLO, the expanded result is given by 

.2 




C{z,M,Hf) =(5(1 - z) ^ -< (5(1 - z) 


TT 


IStt^ Cl ^ 


+ 2CaP[{z) 


+ 


a.. 


TT 


CISa{z) + CATpUfSfiz) + 5i(z) 


where the auxiliary function P^{z) is defined as 


Pniz) = 


l-Z 


In" 




Besides, the two-loop coefficients Si{z) are defined as 


Sa{z) =(5(1 - z) 


Ik 

2 


144 


1 

12 


Lm — 


55 C 3 87 l 7 r 


72 




288 




607 ' 
864 


Sf{z) =(5(1 - z) 




1 


36 6 


5(^(3 657r^ 


'zzTr + X Lm + -rz — 


18 216 


41 


+ Ip'{z)-I^Piiz) 




S'i(z) =(5(1 - z) 


(3.1) 


(3.2) 


+ Psk) - ^P2iz) + (^ - A'(^) + + ^ + 


23 ,,, 21 r 1 /r ^ 21 5l 82 

- (2ci -|- 2(5i -|- -KM Lm H - ((5i + 77r^) ci H - 1 - 

96 ; M T ^ ; 1 ^ 32 ^10 48 


499(^3 


297r^ 2177r2 7451 

+ wwrr + + 


128 


192 


864 



3 37r^ 

-h 

2 (ci + ^i) + — 


P[{k 


(3.3) 


where the notation Lm is defined as Lm = ^T^i^zz/kf)- Tl^® NLO and NNLO results 
obtained this way include contributions from the complete one- and two-loop virtual cor¬ 
rections ^ and soft gluon real corrections. All the scale dependent parts are process inde¬ 
pendent and have been incorporated in various anomalous dimensions. The unknown scale 
independent parts are represented by 5i^2- 


4 Numerical results 

In numerical calculation, we take the SM input 

rufi = 125.7 GeV, mt = 173.2 GeV, = 4.89 GeV, 

Mz = 91.1876 GeV, Tz = 2.4952 GeV, Mw = 80.398 GeV, 

Gf = 1.16639 X 10"® GeV^ a{Mz) = 1/132.338. (4.1) 

^This means that the actual diagrams at NLO and NNLO are of two and three loops if the LO is already 
of one loop. 


































M [GeV] 

41 


Figure 3. The cross section and A"-factor for the pure signal processes at the 8 TeV LHC. 


The factorization and renormalization scales are set to be Mzz- When discussing the scale 
uncertainties, we vary them from M^^/2 to 2Mzz- We use the MSTW2008LO, NLO, 
NNLO PDF sets [41] and associated strong coupling constant to calculate the LO, NLO 
and NNLO results, respectively. We are interested in the results at the 8 TeV and 13 TeV 
LHC. In our calculation, we also consider the decay of the Z bosons, e.g., ZZ e~^e~ , 

so the invariant mass of the four final leptons M^i = Mzz- 

The unknown NLO and NNLO results for the non-logarithmic parts in the hard func¬ 
tion of the interference process, i.e., and 62 , are estimated by using the analogous 
results in the pure signal process and varying by a factor of ^ 1 ^ 2 - More specifically, the 
non-logarithmic parts in the hard function of the interference process are approximated as 

<5i = 6 ^ Y 

<^2 = '^2 C2. (4.3) 

The default value of ^ 1,2 is chosen to be 1. We vary ,^ 1^2 from 0 to 2 to estimate the 
uncertainties coming from those uncalculated non-logarithmic parts in higher order virtual 
corrections. 

Before we present the resummed result for the signal, we first estimate that to what 
extent the resummed result is valid. Since the pure signal process has the same initial and 
final state as the inference process, we take it to illustrate this point. The LO, NLO, and 
NNLO singular contributions, obtained according to eq.(3.1) in the above section, along 
with the exact results, calculated by FehiPro ^ [42, 43], are shown in figure 3. Their ratios 
are shown in figure 4. We first notice that the exact NLO and NNLO iL-factors are very 
significant, and almost constants over a large region of Mzz > ‘^Mz- Then it is evident 
that the contribution of the singular terms dominates the higher order corrections. The 
singular contribution almost coincides with the exact results for Mzz > ‘2-mt at both NLO 

®Both the claimed exact NLO and NNLO results contains exact top quark mass dependence only up to 
NLO. 
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Figure 4. The fraction of the singular terms to the exact cross sections. 




Figure 5. The cross section and uncertainties for the interference processes at the 8 TeV LHC. 
In the left plot, the LO scale uncertainties are about ±20% ^ ±30%, not shown in the plot. The 
uncertainty in the right plot refers to the uncalculated multi-loop amplitudes of the background. 


and NNLO; the difference is below 2% at NNLO. In the smaller Mzz region, the singular 
contributions are a little less than the exact results. But the difference is less than 4% and 
10% at NLO and NNLO, respectively, for Mzz > 220 GeV, which is the off-shell region 
defined in experiments [14, 15]. Moreover, the scale uncertainties of the singular terms lie 
in or highly overlap with those of the exact results for Mzz > 220 GeV, as shown in the 
right plot of figure 3. From this comparison, it is reasonable to use only the singular terms 
to predict the unknown higher order effects. 

We now provide such a theoretical prediction with higher order effects for the interfer¬ 
ence processes between the diagrams 1(a) and 1(b). 
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Figure 6. Same as figure 5, but at the 13 TeV LHC. 

As shown in figure 5, the contribution of the interference process is negative even 
overwhelming the positive pure signal, and has a sharp peak valley around Mzz = 200 
GeV. There is also a small valley at Mzz = 380 GeV, which becomes more significant at 
higher orders. In contrast, the two peaks in the pure signal are of almost the same height. 
We have checked that the shape of the differential cross section actually depends on the 
choice of the scale and kinematical cut. If a fixed scale, e.g., mh is used, then all the 
differential cross section would increase and the cross section in the larger Mzz region gets 
more significant improvement compared to the case of dynamical scales. However, since 
the perturbative expansion of the cross section calculated by a dynamical scale converges 
better, we choose the dynamical scale in our calculation. 

The results at the 13 TeV LHG are shown in figure 6. The shapes of the differential 
cross sections are almost the same as at the 8 TeV LHG, as shown in figure 5. 

The iL-factors for the interference process at the 8 and 13 TeV LHG are shown in 
figure 7. They increase from Mzz = 130 GeV to Mzz = 350 GeV « 2mt, where the 
interference contributions are most significant. Then they decrease with the increasing of 
Mzz, and nearly unchanged for Mzz > 500 GeV. The dependences of the iL-factor on 
the invariant mass at the 13 TeV LHG are similar, but the values are a little smaller. The 
overall NNLO iL-factor is in the range of 2.05 — 2.45 (1.85 — 2.25) at the 8 (13) TeV LHG. 
Here we point out that the ratio of ATnnll'/-^nlo (about 1.3) is similar to that in the Higgs 
pair production we have studied earlier [29], as expected. We also observe from hgure 7 
that ATnnll'/-^NNLO is about 1.1 at both the 8 and 13 TeV LHG. This means that the soft 
gluon resummation is important in providing more accurate theoretical predictions. 

“^We have used the notations NLOsing and NNLOaing to denote the contributions from singular terms at 
NLO and NNLO, respectively. We keep the notations in figures and tables, but neglect the subscript in the 
text for simplicity. 
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Figure 7. The X-factor for the interference process at the 8 and 13 TeV LHC. 
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Figure 8. Comparison of the itT-factor in the pure signal and interference processes at the 8 TeV 
LHC. 

Figure 8 shows that the it'-factors in the pure signal and interference processes are the 
same, which is in agreement with the statements for gg ^ H (600 GeV) —>• WW [3] and for 
gg ^ ZZ without an intermediate resonance [24]. 

Next, we discuss the theoretical uncertainties in the results about the interference pro¬ 
cess, which have been shown in the bottom plots in figures 5 and 6. The scale uncertainties 
at LO, NLO and NNLO are obtained by varying = 1^/ in the range of [M^ 2 / 2 , 2Mzz]- 
For the results at NNLL^, we first vary gf, gs and gh independently by a factor of 2, and 
then combine the errors in quadrature. The LO scale uncertainties are ±20% ~ ±30%, 
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Table 1. The cross section for the interference process in the region 220 GeV < Mzz < 1000 GeV. 
Besides, we also list the scale, parameter 81^2 and PDF+a^ uncertainties. 


^ = 8 TeV 


Ascale[%] 

A«[%] 

ApDF-l-Os [%] 

LO [fb] 

-0.056 

±33.7 

-23.6 

0 

0 

±2.5 

-2.6 

NLOsing[fb] 

-0.106 

±14.5 

-13.6 

±6.7 

-6.7 

±3.8 

-3.4 

NNLOsing[fb] 

-0.129 

±2.7 

-5.5 

±10.3 

-10.3 

±4.0 

-3.7 

NNLL'„/„_,2[fb] 

-0.110 

±10.6 

-6.6 

±7.3 

-7.3 

±3.8 

-3.5 

NNLL'[fb] 

-0.140 

±7.4 

-6.4 

±9.3 

-9.3 

±4.3 

-3.9 


Table 2. Same as table 1, but at 13 TeV LHC. 


^ = 13 TeV 


Ascale[%] 

A«[%] 

ApDF-l-Os [%] 

LO [fb] 

-0.189 

±29.4 

-21.4 

0 

0 

±1.9 

- 2.2 

NLOsing[fb] 

-0.339 

± 11.8 

-11.5 

±6.7 

-6.7 

±3.2 

- 2.6 

NNLOsing[fb] 

-0.407 

± 1.2 

-3.9 

± 10.2 

- 10.2 

±3.4 

-3.0 

NNLL'„/„_,2[fb] 

-0.340 

±10.3 

-7.6 

±6.7 

-7.1 

±3.2 

- 2.8 

NNLL'[fb] 

-0.432 

±7.8 

- 6.8 

±8.9 

-8.9 

±3.7 

-3.2 


which are so large that we do not show them out in the figures. At NLO, NNLO and NNLL', 
they are about ±15%, ±5% and ± 6 %, respectively. Therefore the scale uncertainties are 
significantly reduced after including higher order QCD corrections. The uncalculated multi¬ 
loop amplitudes of the background are evaluated by changing the parameter ^ from 0 to 
2, as described at the beginning of this section. The numerical results are also shown in 
figures 5 and 6 . The associated theoretical uncertainties are about 5% — 10%, depending 
on Mzz, and slightly increase from NLO to NNLO and NNLL'. 

In tables 1 and 2, we list the PDFTa^ uncertainties for the interference process in the 
region 220 GeV < Mzz < 1000 GeV. They are at most about 4%, much less than the other 
uncertainties. If all the theoretical uncertainties are added in quadrature, the uncertainties 
at NNLL^ are about ±12%. The results at the 13 TeV LHG are similar. For comparison, 
we also list the resummation results without tt^ enhanced terms by setting the hard scale 

= M'^ 2 : ill 60 - (2-9), denoted by NNLL'„/o_.„. 2 . We find that this kind of resummed 
cross section is close to the NLO singular terms, much less than the NNLO singular and 
NNLL' results, which means that the vr^-enhanced terms make the main contributions to 
the NNLL' results. This feature is in agreement with that in single Higgs production [30] 
and the double Higgs production [29] at the LHG. The scale uncertainty at NNLL'^/o _^2 
is reduced compared to the NLO singular terms, but a little larger than total resummed 
result. The other theoretical uncertainties at NNLL'.^,,/o _^2 are nearly the same as the NLO 
singular terms. 
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Figure 9. The ratio of the contribution from the interference between the diagrams 1(a) and 2(c) 
to the LO cross section. 

Now we evaluate a special contribution in gg{-^ h*) —)■ ZZ but not in gg{-^ H) —>■ 
WW, i.e., the interference between the diagrams 1(a) and 2(c), which would appear at 
the NLO corrections. Furry’s theorem states that only the axial vector part of the Z 
boson coupling can contribute. Since the coefficient of the axial current is proportional 
to the weak charge T ^„3 of the SU{2)l gauge group, the contribution is proportional to 
the difference of of the quarks in a SU{2)l doublet. Thus only the third generation 
quarks, massive t- and b-quarks, generate non-vanishing result, which is both infrared and 
ultraviolet safe. As shown in figure 9, they are so small that we can neglect them when 
considering higher order corrections. 

5 Conclusion 

We have studied the high order QCD effects, in particular the soft gluon resummation, in 
the signal-background interference process of gg{-^ h*) —)• ZZ. This process can be used 
to constrain the total width of the Higgs boson and provide a special way to measure the 
couplings between the Higgs boson and other particles in SM. The previous theoretical 
prediction for the interference process is only at the LO because of technical difficulty in 
calculating massive two-loop integrals. In this work we show approximate NLO and NNLO 
cross sections obtained from the resummation formalism in SCET, and also present the 
result with the soft gluon effect resummed to all orders in Comparing the approximate 
results with the exact ones for the pure signal process at NLO and NNLO, we observe that 
approximate results almost reproduce the exact ones in the off-shell region, and find that the 
high order QCD effects for the interference process are very sizable. For the interference 
process, the approximate NNLO AT-factor is in the range of 2.05 — 2.45 (1.85 — 2.25), 
depending on Mzz, at the 8 (13) TeV LHC. Besides, the soft gluon resummation can 
increase the approximate NNLO result by about 10% at both the 8 TeV and 13 TeV LHC. 
At the same time, the corrections from the soft gluon resummation are similar to that in 


- 14 - 



the Higgs pair production. We also find that the approximate ii"-factors in the interference 
and the pure signal processes are the same, which is in agreement with the statements in 
previous literatures [3, 24]. 

Moreover, we discuss the theoretical uncertainties in the results. The scale uncertain¬ 
ties are significantly reduced after including higher order QCD corrections. The uncer¬ 
tainties from uncalculated multi-loop amplitudes of the background slightly increase from 
NLO to NNLO and NNLLk The PDF-t-a^ uncertainties are rather small compared with 
the others. Combined all together, the theoretical predictions at NNLL' suffer from roughly 
0(10%) uncertainties. 

Our study can be easily extended to other processes, such as the pure signal process 
gg ^ h* ^ VV and pure background process gg —)• VV iV = VF, Z) via quark box-loops. 
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